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Abstract 

I I We present a new study by the STAR Collaboration of background 

^ fluctuations in jet reconstruction in heavy ion collisions. 

' I Jets are fundamental objects in QCD, and are ubiquitous in high energy 

O colhsions of aU kinds. In heavy ion collisions wc utilize jets to probe the hot 

2 medium. While the rate of initial hard (high Q^) scattering is not altered by the 

I I presence of medium, the interaction of the scattered partons with the medium 

induces changes to the internal jet structure, and may alter the jet direction. It 
is the modification of internal jet structure and of inter-jet correlations that are 
the observables of jet quenching [1]. 

Initial studies of jet quenching at RHIC, utilizing high px hadrons and their 
correlations as jet surrogates, have found marked suppression effects that estab- 
CN lish jet quenching experimentally [T] . Extensive effort is underway to understand 

these measurements theoretically T, ^ and to use them to quantify medium 
T-H transport properties 0^ . While such hadronic measurements have been valuable 

to establish the phenomenon of jet quenching and provide initial quantitative 
" , I constraints on medium properties, their scope is limited by the biases inherent 

^ in leading-hadron observables [1], i.e. the observed high px hadron popula- 

tion is dominated by remnants of those jets that have interacted least with the 
medium. Full study of the dynamics of jet quenching requires unbiased mea- 
surement of jets that have undergone significant interaction in matter, either via 
jets recoiling from a direct photon or Z boson, or via unbiased jet reconstruction 
that recovers the energy "fiow" of the scattered parton, independent of how it 
fragments. 

The measurement of jets in high energy heavy ion collisions is more complex 
than in elementary collisions, due to the large fluctuations in the underlying 
background. The first measurement required in such a program is the inclusive 
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jet cross section, which quantifies the degree to which unbiased jet reconstruction 
has been achieved in practice. Prehminary STAR measurements of the inclusive 
jet cross section at mid-rapidity in central Au+Au collisions at ^/snn— 200 GeV 
have been reported previously 0, based on the kr and ant i-Zcp [Hill] algorithms 
and the Fast Jet background correction scheme |H1 STAR measures jets in 
its central barrel (|?7| < 1.0) utilizing charged particle tracking and electromag- 
netic calorimetry [TUj. Underlying event fluctuations distort the inclusive jet 
distribution significantly, and the corresponding correction is by far the largest 
source of systematic uncertainty in this measurement. 

The Fast Jet-based background correction to the inclusive cross section mea- 
surement proceeds in two steps: 

1. Event-wise background estimate: reconstruct all jets in the experimental 
acceptance using the fcy algorithm with R = 0.4 [71 [5] and determine the 
median energy density p = median{p^rp / A^} , where and A'- are the pr 
and area of the ith jet in the event, 1 < ? < Njef Corrected jet pr is then 
p-hcorr = _ p . This accounts on average for the underlying event 
background, but not residual fluctuations. 

2. Deconvolution of fluctuations: estimate inclusive spectrum distortion due 
to fluctuations; express in terms of bin migration and correct by regular- 
ized matrix inversion. 

In [S], the spectrum distortion in step (2) was estimated by generating 
PYTHIA jets with a realistic pr spectrum for p + p collisions at y/s=200 GeV, 
injecting (or "embedding") them into real STAR central Au + Au events, and 
reconstructing the hybrid events with the same procedure as used for the data 
analysis. The resulting spectrum distortion was parameterized as the convo- 
lution of a Gaussian function with a = 6.8 GeV and systematic uncertainty 
±1 GeV. While this is a well-controlled procedure, it does not address two im- 
portant issues: (i) sensitivity of the correction to the fragmentation pattern 
of the signal jets, of particular concern since fragmentation of quenched jets 
is not modeled by PYTHIA, and (ii) non-Gaussian response, especially a non- 
Gaussian tail towards positive fluctuations that can distort the spectrum over 
a wide pt range. 

In these proceedings we present a new study that addresses these issues. 
This study is likewise based on embedding of a simulated jet population into 
real STAR central Au + Au events, but with a wider spectrum of jet models than 
the previous study: single hadrons, PYTHIA jets, and Q-PYTHIA jets to model 
quenching effects [TTj . Embedding is carried out at the level of charged particle 
tracks and calorimeter energy deposition, without accounting at present for 
instrumental effects in order to isolate the irresolution due solely to background 
fluctuations. For the background population we use 8M STAR Au + Au central 
(0-10%) events, from an event sample recorded with a minimum bias trigger in 
the 2007 RHIC run. 

We embed an object of known pt^Pt"'"^'^ and apply jet reconstruction to 
the hybrid event, utilizing the FastJet anti-fc^ algorithm with R = 0.4. We 
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find all reconstructed jets containing embedded charged tracks and calorimeter 
towers, and identify the matched jet as that reconstructed jet containing over 
50% of the total pt of the embedded jet. The matching efficiency is found to 
be high (over 90% at low p^, over 95% for pT>lQ GeV). For embedded p'^^^'^ 
and matched jet p^'^° , we quantify the response of the hybrid system to the 
embedded jet via: 

Xrr, o^'''^co Areco embed /■•i \ 

where A'^'^''° is the area of the matched reconstructed jet and p is determined 
prior to the embedding step. This definition is identical to Eq. (1) in [5|. The 
normalized distribution of 6pT is the probability distribution to find jet energy 
(after event-wise background correction) p^f^"^ = p*™'^ + SpT- If there were no 
background fluctuations, 5pT would be a delta function at zero. 




Figure 1: Single particle embedding: distribution of ^^^(Eq. [T]) for a "jet" 
corresponding to a single pion with pT=30 GeV embedded into central Au + Au 
events. The same dpT distribution is shown in both panels. Left panel shows fit 
with Gaussian function, right panel shows fit with Gamma function modeling 
independent particle emission [IJj. See text for details. 

We first consider the embedding of a single hard pion. While this is an un- 
likely fragmentation pattern for a jet, its simplicity aids in interpretation and 
we find that single particle embedding captures many of the features of the sys- 
tem response. We will later compare the single particle response with that of 
more physical jet fragmentation patterns, in order to isolate the specific role of 
fragmentation. Fig. [T] shows the normalized distribution of 6pT for a single pion 
with pt=30 GeV, rj — 0.0 and (f> = (i.e. far from STAR acceptance bound- 
aries), embedded successively in 8M different STAR central collision events. It 
is the ensemble-averaged response to this particular choice of jet probe. The 
probability to reconstruct a 30 GeV single-particle jet with greater than 20 GeV 
in excess energy (after event-wise background correction) is ~ 10"'^. In order 
to characterize the distribution we fit with two analytic functions: 

• Fig. [T] left panel, shows a Gaussian function fit to the region of SpT less 
than the median of the distribution. The fit gives a ccntroid -1.27 GeV 
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and width 6.17 GeV. The region Spr < is well modeled as a Gaussian 
but the region Spx > has significant excess relative to the fit, which 
can strongly influence the inclusive spectrum measurement. A Gaussian 
function is not an accurate model of the distribution of fluctuations [5] • 

• Fig. [T] right panel, shows a fit to the region Spx < of an analytic func- 
tion describing the transverse energy distribution into finite acceptance of 
multiple independent particle emission sources, each emitting an exponen- 
tially falling spectrum [T^]. The spectrum has fixed (pr) — 0.5 GeV/c, 
corresponding roughly to RHIC conditions, and the fit finds about 370 in- 
dependent "sources". The fit function is non-Gaussian at large Spx- This 
simple model of uncorrelated particle emission can account for the bulk of 
the irresolution due to background fluctuations. 




Figure 2: Jet area vs. Spx for embedded single pions with pT=0.l (upper left), 
5.0 (upper right), 10.0 (lower left) and 30.0 (lower right) GeV. 

Background fluctations are generally expected to be correlated with iet area; 
for background due to uncorrelated emission, ajiuct ^ VA [311]. Fig. shows 
the distribution of measured jet area [7l[8] vs. Spr, for embedded single pion jets 
with a broad range of pr- 0.1, 5, 10 and 30 GeV. At large px, there is a narrow 
area distribution around the expected value ttR"^ ~ 0.5 for R = 0.4, with little 
dependence of the 6pT distribution on pT- At lower px there is a significant 
population of lower area jets, due to the response of the anti-fcT algorithm in 
dense background for jets lacking an energetic core. The Spx distribution is 
narrower for smaller jet area, as seen in Fig. [3j left panel. The deconvolution of 
background fluctuations from the inclusive jet spectrum must accurately take 
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into account the pT-dependence of the jet area distribution, and consequently 
the pT-dependence of fluctuations (i.e. of Spx)', this is the subject of forthcoming 
work. 




Figure 3: SpT distribution systematics. Left: for different jet areas, for probe 
with pt=0.1 GeV. Right: overlay of Spx for single pions, PYTHIA and Q- 
PYTHIA jets, of varying px- 

Finally, we investigate the dependence of Spx on fragmentation pattern and 
Pt- Fig. [Sj right panel, shows the overlay of multiple Spr distributions for 
single particle jets and for jets generated by PYTHIA and Q-PYTHIA {q=5 
GeVVfm). For both PYTHIA and Q-PYTHIA, 30 different jets within each 
of the intervals 1 < pt < 5 GeV and pr > 30 GeV were embedded. In order 
to compare their shapes directly, the distributions were aligned horizontally by 
fitting a Gaussian function to dpT< and aligning the centroids by shifting 
relative to one reference distribution. The shifts are shown in the insert and 
are typically smaller in magnitude than 1 GeV. The overlay shows that the Spx 
distribution is to a large extent universal, within a factor '-^ 2 at SpT=50 GeV. 
Two energetic jets from the Q-PYTHIA population do not follow the bulk of 
the distribution for Spx < — 10 GeV, though this will have little impact on the 
smearing of the inclusive spectrum. 

Insensitivity of the Spx distribution to fragmentation pattern is crucial for 
well-controlled analysis of quenched jets, whose fragmentation is a priori un- 
known. Fig. [3j right panel, indicates that reconstruction in a dense background 
environment using the anti-fc^ algorithm [B] is indeed robust against variations 
in fragmentation. Further quantification of this observation, and its application 
to deconvolution of the measured inclusive jet spectrum in central Au + Au 
collisions, is in progress. 
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